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1. Introduction 

The Bose-Hubbard Model with nearest-neighbour hopping is defined on 
a cubic lattice Z rf of dimension d by the Hamiltonian 

H=- (a* x -a* y )(a x -a y ) + X^2n x (n x -l), (1.1) 

x,y&Z d : \x-y\=l x 

where the operators a* and a x are creation and annihilation operators sat- 
isfying the usual bosonic commutation relations 

\(l Xl dy\ — &x,y 

The operators n x = a* x a x are the local number operators. 

In pQ a long-range-hopping version of this model was analysed. It is given 
by the Hamiltonian 

x v 

Hv= 2V S ^ ~ a l^ ax ~ a v> + A Yl n ^ nx ~ l ) ( L2 ) 

x,y=l x 

on a complete graph of V sites. In particular, the following variational 
expression for the pressure was derived: 



p(/3,fj,,X) = 

sup | — r 2 + — InTr exp + A — l)n — An 2 + r(a + a*)) j. 

(1.3) 

The derivation made use of the so-called approximating Hamiltonian method 
introduced by Bogoliubov Jr. [3j and made rigorous by Zagrebnov et al. [3], 
see also [5]. Here we present a new derivation of this formula using the C*- 
algebraic method of Petz, Raggio and Verbeure [2j, which is inspired in part 
by Varadhan's theorem in probabilistic large deviation theory, and in part by 
work of Fannes, Spohn and Verbeure [B]. Our C*-algebraic approach corre- 
sponds to variational expressions for thermodynamic functionals of classical 
Gibbs measures [7], and in particular we are using the variational expres- 
sion for the relative entropy (see [8] for classical probability theory). As a 
result we obtain a variational expression for the pressure where we optimize 
over states of the infinite system. Due to the symmetry (lack of geometry) 
on the complete graph this expression can be simplified to (|1.3p using the 
well-known St0rmer theorem. Hence, our results here are a quantum anal- 
ogy of recent results for classical statistical mechanical models on complete 
graphs using exchangeability [Qj . The analysis in [2] concerns quantum spin 
models and the extension to the Bose-Hubbard model requires a number of 
technical considerations. Many of these can be found in the book by Ohya 
and Petz |10j . but we present some of the proofs here nonetheless in order 
to make this paper more self-contained. 
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Section [2] contains the proof of our variational formula (|1.3p . In Section [3] 
we briefly present some of the features of the model again, mainly in order to 
correct some minor but irritating errors in the analysis of [l] . In Section HI 
we show that some of the C*-algebraic formalism extends to the nearest- 
neighbour hopping model, resulting in a variational formula for the pressure 
analogous to the well-known formula for spin models. This formula does 
not appear to have been written down before. As the variation is over the 
set of all translation-invariant states on the lattice, it is difficult to analyse, 
however, just as in the case of spin models. 

2. C*- ALGEBRAIC DERIVATION 

We follow the technique of Petz, Raggio and Verbeure [2j . The Hamilton- 
ian is invariant under the permutation group. We write it as 

x=l y=l x=l 

where 

h x = n x + Xn x (n x - 1). 
For finite V we thus assume that a CCR algebra Ay is given generated by 
creation and annihilation operators a* and a x (x G {1,...,V}) and with 
standard representation on the Fock space Ty . 

Next we define a reference state ojy on the Fock space Ty as the product 
state ojy = (2)^=1 ^x, where u> x has the density matrix 

_ exp[/3Qre x - h x )} 
Tr exp\J3(fm x - h x )} ' 

Let 

oo 

A = (J Ay (2.2) 
v=\ 

be the (quasi-local) CCR algebra for the complete lattice N, where the clo- 
sure is taken with respect to the norm topology. 

Lemma 2.1. Suppose that 4> is a regular permutation-invariant state on A 
such that (j)( n x) < +oo for all x 6 {1, . . . , V}. Then 

1 V 1 V 

x=l y=l 

Proof. Notice that we can define (f>(n x ) as the supremum 

(j){n x ) = sup (f>(pffln x ), 

N>1 
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(x) 

where is the projection on the subspace of T x with n x < N. Ob- 
viously, this is independent of x by permutation invariance. The formula 
^((f T^c=i a x) {y Yly=i a y)) should be interpreted in a similar way: 

1 V 1 V 1 V 

x=l y=l — x,y=l 



where 



x<=V 

We now write 

v , v 



((^E<)(^E% 

x=i j/=i 

1 V 



x=l x=£y 

The first term is clearly bounded by y(/>(n x ) and hence tends to zero. By 
permutation invariance, the second term equals 

^ ~~ 1 < '^v „*„ r>V\ 



V 



{P%ata 2 pV) 



We conclude by proving that the limit (first N — ► oo and subsequently 
F — > oo) of this expression exists and equals 



<t>(a\a 2 ) = lim <j>{P^ a\P^ P^ a 2 P^). (2.3) 

iV— >oo 

To this end we write 

= ((P& - P^K^J + (P^a 2 (P^ - p£j) 

and treat each term separately. Both terms are similar; we consider only 
the first. We have 

l0((^-O«i«2P^)i 

< [H( p m -P&W'iiPl ~ p nM' 2 [H p nA^pI)} 1/2 ■ 

The second factor is obviously bounded by ^{n^) 1 ^ 2 ■ In the first factor we 
can write 

+ ... + <«...®<-"® ( pi;; ) -p«r ) )- 
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With the observation that — P^ 2 commutes with a\a\ we get 

< - P$)a* iai (pW - pW)) + (V - - Pg), 

and therefore the limit of P^a%a y P^ exists. Similarly, one proves that 

^P^ala 2 Pl) - ^pUatpUpVazPW] 

for fixed V. Taking the limit N — > oo and subsequently V — ► oo the result 
follows. □ 

Remark 2.2. Since is regular, its restriction <j>v to each Ay is regular, and 
the number operators n x are well-defined. Moreover, the corresponding GNS 
representation is equivalent with the Fock representation by Von Neumann's 
theorem [11]. In particular, cfty is normal for all V, i.e. <j> is locally normal. 
Thus 4>y has a density matrix p^. 

Lemma 12 . 1 1 gives the mean energy for our model. Next we are concerned 
with the relative entropy with respect to our reference state ojy. It is well- 
known that the relative entropy 

S((f)v || uy) = Tr [p0 v (lnp 0v - lnp^)] 

is convex and superadditive [10j . A precise definition of S(<p \\ u) for states 
on a Von Neumann algebra was given by Araki [12, 13J. Given the stan- 
dard representation n of £>(Py) (the GNS representation with respect to the 
tracial state) one has 

S(0y \\0Jy) = -($y| IllAn v ,* v \$v), (2-4) 

where 

1/2 1/2 

is the relative modular operator. Here <3>y = pi and £ly = pj v respectively, 
and hence 4>y(A) = ($y \ir(A)$v) and ujy(A) = (Qy, 7r(A)Qy). If Pn,$ is 
the corresponding resolution of the identity, then we can write 

S{4> v \\lo v ) = - [ ln\(<S> v \E Qv ^ v (d\) |$y) 
J 

/oo 
lnA(*y|Pn F)$F (dA)|$y). (2.5) 

Here we have separated the two integration domains to indicate that the 
second integral is always convergent, whereas the first should be interpreted 

as 



-inf / lnA($y|Pn v ,$ v (dA) |$y) 
5 >o Js 



and determines whether S(</>y || uy) is finite or infinite. Note that 
S{4>v || wy) = +oo if $y j£ supp (cjy). 
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An equivalent definition was introduced by Uhlmann 

S{cf>y \\uv) = -Hm~ (||A<£ ^ v f - l) . (2.6) 

This is easily seen to be equivalent using Lebesgue's monotone convergence 
theorem. We now prove the following variational formula for the relative 
entropy (compare [SJ Lemma 3.2.13] for classical probability theory): 

Theorem 2.3. For any state <j>y on B{J-y), 

S{fo\\u v )= sup |/30y (A) -In Tre^-£-i^ +A )) 

AeB{F v ): A*=A 1 J 

+ V In TV e l3 ( fJ ' ni ~ hl \ 

Proof. We first prove that S{(f)y \\ ujy) is greater than the right-hand side. 
Let A € B(Ty) be self-adjoint. Then we can define the perturbed state tpy 
by 

e /3(Mf-ELi ft»+A) 
^ v ~ Tre^^-ELi^+^)' 

We write 

Z A , y = T re^-^i^ +A ) 

and Zy = Tr qP^ j ^'~^'x=x h * in the following. We employ a change of state 
(measure) method [8] with respect to the reference state coy and the per- 
turbed reference state rpy. The non-commutativity of our random variables 
is an additional difficulty. But the Du Hamel formula gives 

(*v|A^ Vj * v |*v) 

= Z^Z v ^v\ (l-/3^ t dr7r(e-^A^)7r(p^)7r / (^)|$ V/ 

(2.7) 

where 

v 

x=l 

and where H v = —fiAf + Ylx=i h-x is the non-interacting part of the Hamil- 
tonian in (|2.ip . Differentiating we get 



A| ($ y |A^ v $v = -lnZ A>y + lnZy-/3($ V '|7r(A)$v') 

d (2-8) 
+ ^Uo+^l A W,*vl # v) 

and hence 

S(<p v || W) = II wv) - PMA) + In - ^ v — . (2.9) 
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The desired inequality now follows from the positivity of the relative en- 
tropy. (Notice that this follows immediately from In A < 1 — A and 

J A($y| Enfy t ^y(dX) |$y) = ($y| Avp Vj $ v |$y) 

= (^^k(%)f'Wv)^ |*v) = Tr(p^) = 1 

by a simple approximation.) 

To prove the converse inequality, first assume that c\ojv < </>y < C2^y 
for constants < ci < C2 < +oo. Then there exists a bounded relative 
Hamiltonian A such that 



c 



Tre /3(MAT-ELi^+^)' 

that is, 



Indeed, it follows easily that Dom(ln/?^ v ) = Dom(ln/9^ v ) and A = In, 
In / o a;v is bounded. The identity (|2.9|) with 0y = -i/v then yields 



S{(f> v \\u v ) = /3^ v (A)-laTie^-^=^ hx+A ^ + VlaTre' 3 ^ ni - hl \ (2.10) 

The general case then follows from the lower semi continuity of the relative 
entropy [10] : 

S(<pv II wy) < liminf 5(0y e || wv). 
y— >oo 

Indeed, first assuming 0y < Xujy we can put cfiv,e = (1 — e)0y + to 
conclude that the theorem holds in this case. In the general case, we use 
the approximation p^ Vt = j^pf^ p e ] • □ 



It is proved in [10] , Corollary 5.21, that 

S{(j) || Wi ® W2) > 5(01 || Wi) + 5(02 || W2). (2.11) 

The relative entropy is therefore super additive: 

5(0yi+y2 || WVi ® Wy 2 ) > 5(0Vi || W^) + 5(0y 2 || uy 2 ). (2-12) 
It follows that the mean entropy 

s(0||w):= lim —5(0y || wy) 
y^oo y 

exists. We now have the following 'level-Ill' variational expression for the 
pressure. Here level-Ill refers to the fact that in the variational formula 
we optimise over states of the infinite system [7J. The set of all regular 
translation-invariant states on A is denoted by S(A), and the set of all reg- 
ular translation-invariant and permutation-invariant states on A by Su(A). 
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Theorem 2.4. 



p(/?,M,A) := Jim -^InTre^^ 
»oo pi/ 



sup (0(a*a 2 ) - \s{(f) || a;)) + \ In Tr< 

i(=ST, : P ) P 



_ • e /3(/i«i-'ii) 

<H 71 3i)< + 00 



i7ere i/ie supremum is taken over all regular translation- and permutation- 
invariant states (ft on A such that 4>(n x ) < +00 for all 

Proof. We denote 

1 V 1 V 

v v=(vT, a x)(vT, a y 



x=i y=i 



and 



so that 



p(f3,/j,,\) = lim \-Py- 



We approximate vy using a cut-off, and call this bounded operator also vy- 
By Theorem[ 



S((ft v \\u v ) > pV(ft v {v v ) -lnTre^-^*=i h * +Vvv) + VlnTre l3 ^ ni - hl) . 
This implies, using Lemma 12. 11 that 

liminf -^-Py > sup \(ft(ala 2 ) - \s((ft \\ lo)\ 



_|_ _ ^rpr e /3(M"l-'il) - 



To prove the converse we seek an approximate maximiser. This is stan- 
dard. We let 

iftv = ipv ® ® • • • 
be the infinite tensor product of states and define 

1 y 

0=1 
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where t~ is the translation over j. This is a permutation- invariant state. We 
estimate the expectation of the energy density: 

1 - 

VvK^) = — ^2i>v(at a k+i) 



V 

k=l 
V-l 



y ^2 Vv(4°fc+i) + y^v(av)^v(ai) 
k=l 



x v v 1 
72YI Yl ^v(a* x a y ) + —i) V {av)^v{ai) 



(2.13) 



y2 /L^i Z-^i K-x-yj < y 

x=l y=l;y^x 
V 



= ^2 ^v(a* x a y ) - — Vv(o*oi) + ^l^y(ai)| 2 . 
x,y=\ 

Lemma 12. II and the Cauchy-Schwarz inequality |?/v(ai)| 2 < i/v(o*oi) then 
imply that 

^(0*02) -ipv(v v )\ ^0 as V ^00. (2.14) 



It is known that the entropy is convex in both arguments [TO]. In partic- 
ular, we have 

s(A<£i + (1 - \)4> 2 || u) < As(0i || u) + (1 - A)s((/> 2 || w). 

On the other hand, by a simple approximation, we have 

S(X<f> iy + (1 - \)(f> 2 ,v II wv) = - A($i,y| In A* Vj n v |$i )V ) 

- (1 - A)($ 2 ,v| lnA^jv |*2,v>, 

where 4>v = ^4>i,v + (1 — X)<p2,v- Using the fact that n(pn v ) and Tr'(p<f> v ) 
respectively tt'(p<j )1 v ),7r / (/)^ 2 v ), commute and the operator monotonicity of 
the inverse and the logarithm, we have 

In A Q < IxlAq v ^ 1v - In A. 

Of course, a similar inequality holds w.r.t. 02 v- Therefore 

S(^<j>i,v + (1 - A)^2,y || wy) > AS^iy || U7) + (1 - X)S((p2,V II Wy) 

+ A In A + (1 - A) ln(l - A). 

In the limit, we find in combination with the convexity above, that the mean 
relative entropy is affine in the first variable: 

s(A0i + (1 - A)0 2 || u) = As(0i || uy) + (1 - A)s(0 2 || w). (2.15) 

Applying this to the state ipy we get 

1 V 

s(ip v || w) = tf^ 5 ^ or fc-iM 
fc=l 



10 



STEFAN ADAMS AND TONY DORLAS 



provided the right-hand side exists. However, by the translation-invariance 
of u, the right-hand side can be written as 
! V 1 V 

T7y^ s (^V °Tfc-ilM = hm — -pr V]5((^y)® n OTk-i\\ w„y) 

k=l k=l 

= hm J-5((^v)®» || (uv)®*) 

n— *oc nV 
= yStyv \\UV)- 

We therefore have 

(Vv(a*a 2 ) - ^ S 0/V II ^)) - (Vv(wy) - -^S(ipv II -> (2.16) 

as V — > oo. On the other hand, by Theorem 12.31 as above (see Eq. (|2.10p ). 

T ie ^-T,Li^+Vv v ) 
S(W || uy) = PVMw) ~ In Trem/mi _ fel) (2 1?) 

= P^v(vv) ~ fiPvPV + lnTre^' mi -H 

□ 

The variational expression for the pressure can now be simplified by de- 
composing the state (f> into an integral of extremal permutation-invariant 
states: 

Theorem 2.5. 

KAjM) = sup {Ka0| 2 - ^SfaWui)} + l-hxTreP^-W, 

cr(ni )< + oc 

(2.18) 

where the supremum is now taken over regular states a of A\ such that 
o~(ni) < +00. 

Proof. The set of permutation invariant states S n (A) is a convex compact 
set in the weak*-topology, and it is metrizable because A is separable. By 
Choquet's theorem [T7], we can therefore decompose an arbitrary state <f> 
into an integral 

Jext(S n (A)) 

over the extremal points of S W (A). Here (j, is a probability measure on 
ext (S 7r (A)). But, by St0rmer's theorem [15], the extremal permutation 
invariant states are the product states 

ipa = o ® a (g> . . . , 

where a is a state of A\. Thus, 

ipc (J.(da). 
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Moreover, since (f>(n x ) < oo, we have that cr(m) < oo for //-almost every a. 
It follows that 

0(a*a2) = y <j(ai)cr(a 1 ) /i(dcr) = / |cr(ai)| 2 //(da). 

For the entropy term we use the following lemma [171 Lemma 9.7]: 

Lemma 2.6. Suppose that X is a compact convex subset of a locally convex 
topological vector space. Let / : X — * M be an affine, lower semi continuous 
function on X, and suppose that fx is a (Radon) probability measure on X , 
xq = J x/j,(dx). Then 

J f{x)n{6x) = /Oo). 



Since the relative entropy (and hence the mean relative entropy) is lower 
semi continuous and affine, the lemma applies and we have 

s{cj)\\uj) = j s^o- \\uj) n(da). 

However, since u is also a product measure, s(ip a \\u) = S(a\\u>i). The 
Theorem now follows. □ 

Remark 2.7. Notice that the subset of regular states is also closed in the 
set of all states on A by the Banach-Steinhaus theorem. Indeed, if 4> a is 
a net of regular states converging to in weak*-topology then (f> a (W(tf)) 
converges uniformly on compact sets t G [a, 6], where W(tf) is the Weyl 
operator for / € C^. 



In the following we write a instead of a\ and n for n\. Our variational 
expression for the pressure can be further reduced to 

Theorem 2.8. 

p(fl,H,\) =sup{|z| 2 -/(z)} + ilnTre /3( ^- /l) , (2.19) 
where the rate function I{z) is given by 

( 1 Q^{^n-h+(va* +va)) ^ 

10) = sup {vz + US-- In Tve ^ n _ h) }. (2.20) 

To prove this, we first need a lemma: 
Lemma 2.9. Denote 

p{v) = 1 In TV ^n-h+^+ua)) _ ( 2 21 ) 

This function is convex and satisfies 

p(y) ~ |z^| 4/3 as \u\ -> 00. (2.22) 
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Proof of Lemma 12.91 Differentiating, we have 
where 

H(v) = —fin + h — va* — ua. 
Differentiating again, we get, as in [I], with v = x + iy, 

q2 . Q Q . 2 

Q^P( V ) = \d^ + dB) ^ = P( a + a* - (a + a*)\a + a* - (a + a*)) H(u) , 

Q2 . Q Q . 2 

d^ P ^ = ~\d^~d&) P ^ = P{ a ~ a * ~(a -a*) \a-a*-{a-a*)) H{v) , 
and 

d 2 , ./ <9 d \f d d \ 



dxdy \dv dv )\dv dv , 

= i/3(a + a* - (a + a*) \ a* - a - (a* - a)) H r u y 

It follows by the Cauchy-Schwarz inequality that the corresponding matrix 
is positive-definite. 

To prove the asymptotic behaviour, we first remark that 

zai* + Pa < 2\v\(n + 1) 1/2 

and hence 

p(iz) < ilnTre^^- h+2 ^l(" +1 ) 1/2 ) 

1 oo 
P n=0 

For the reverse inequality we use one half of the Berezin-Lieb bounds 
Tr e -HW > [ dzdE c -(z\BM\z)_ 



7T 

Here \z) stands for the coherent state 



.12 ^ z n (a*) n . . 

n=0 

Since 

(z| a |z) = z, (z| a* |z) = f, 
(z\ n \z) = \z\ 2 , 

and 

(z\(a*) 2 a 2 \z) = \z\\ 

we get 



□ 
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Proof of Theorem 12.81 Take A to be an approximation of va* + va in 
Theorem 12.31 We get, writing 



p(y) = IlnTre^"-^^*^), 

P 



the inequality 



|^(a)| 2 -~S(0y || oj v ) < |0(a)| 2 -0(zaz* + Va) + p(u) - p(0) 
and since v is arbitrary, 

|0(a)| 2 -~S(0 || u)< |0(a)| 2 -/(0(a)). 

Conversely, suppose that z is a maximiser for sup{|z| 2 — I(z)}. For any 
v G C, define the state 0^ by 

g/3(/wi— /i+(i^a*+Pa)) 
~~ r JY e /3(^tn-h+(^a*+Pa)) 

and choose z/ such that (j) v {a) = z. It follows from the above lemma that 
such v exists. Then 

sup{|0(a)| 2 -is(0|M} > |0„(a)[ 2 -S(&,||a;) 

= |z| 2 - 4> v (va* + Va) + p(v) - p(0) 
= \z\ 2 -I(z). 

□ 

We finally rewrite the expression in the form (|1.3|) . With a gauge trans- 
formation it is easy to see that I{z) only depends on |z| and we have 

p(P,fi,X) = sup{x 2 - I(x)\ +p(0) (2.23) 

x>0 

and 

I(x) = sup{2rx -p{r)} +p(0). (2.24) 

r>0 

Now let r > be given, and suppose that x max is a maximiser of the expres- 
sion (12T2^ . Then 

p(/3, /i, A) = x 2 ^ - J(x max ) + p(0) < xl mx - 2rx max + p{r) 

and choosing r = x max , 

p{P, V, A) < -a4 ax + p(x max ) < sup{-r 2 + p(r)}. 

r>0 



On the other hand, if rg is a maximiser of the right-hand side, then 

r °=2^ (r)| -o 



and hence 

/(r ) =2r 2 -p(r )+p(0). 
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Inserting, we get 

sup{-r 2 +p(r)} = -rl+p{r Q ) = - J(r ) + p(0) < p(/3, fi, A). 

r>0 

Hence, we have shown the formula (jl.3p for the infinite range Bose-Hubbard 
model ([23J) . 

3. Analysis of the phase diagram 

The phase diagram of the model was analysed in pQ. The same model, 
but with disorder, was analysed in [16] where it was found that the disorder 
gives rise to new phenomena. Unfortunately, [1] contains a few errors, which 
we wish to correct here. First of all, the critical values of lambda are not 
given by (2.14) of pQ, but instead 

X Ctk = 2k + l. (3.1) 

This was already remarked in [16], see Remark 4.1. Indeed, although a gap 
exists for A > A^ given by (2.14) in pQ, the limiting value of fi((3,X) lies in 
this gap only if A > A Cj fc. 

There was also a mistake in the program to compute the p — V diagrams 
of Fig. 5 and 6 in pQ, as well as the condensation fractions of Fig. 7 and 8. 
We include corrected graphs below: 




The P- V diagram for A = 5 . 
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The condensation fraction as a function of the density for A = 5. 



BOSE-HUBBARD MODEL 



17 



4. The nearest-neighbour hopping model 

As in the case of quantum spin models, there is also a variational formula 
for the pressure of the translation-invariant Bose-Hubbard model, analogous 
to Theorem 12.41 

Theorem 4.1. The pressure of the nearest-neighbour hopping Bose-Hubbard 
model is given by 

d 

p((3, /X, A) = SUp <K a O a c„ + a e„«0) - -aS{(j> \\ Uj) \ 

4>:<p(n x )<+oc L u=1 P } (4 -Q 

+ iln r Tre /3( ^- /lo) , 

P 

where the supremum is over all regular translation-invariant states 4> on 
the CCR algebra such that 4>{n x ) < +oo, and to is the product state uj = 

n - J_ p j3(iJ,n x -h x ) 

^0 

and 

h x = dh x — Xh x (h x — 1). 

The derivation of this formula is completely analogous to that of Theo- 
rem 12.41 The main difference is that the infinite-volume limit now has to 
taken in the sense of Van Hove. For the case of spin models, see for example 
[19j or [20]. This variational formula does not seem to have been written 
down before, though it has to be said that it is not clear how useful this 
formula is. The analogous formula for spin models has so far not been very 
useful for analysing the phase diagram. One possible application is perhaps 
the cluster variation approximation, see [21], [22J, [23], [24] . 
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